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Abstract 

For the first time a rigorous quantum treatment of the Landau- 
Pomeranchuk-Migdal effect in QED and QCD is given. The rate of photon 
(gluon) radiation by an electron (quark) in medium is expressed through 
the Green's function of a two-dimensional Schrodinger equation with an 
imaginary potential. In QED this potential is proportional to the dipole 
cross section for scattering of e + e~ pair off an atom, in QCD it is propor- 
tional to the cross section of interaction with color centre of the color singlet 
quark-antiquark-gluon system. 
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The effect of multiple scattering on bremsstrahlung (the Landau-Pomeranchuk-Migdal 
(LPM) effect in QED and QCD has recently attracted much attention [|,g. How- 

ever, a rigorous treatment of the LPM effect for an arbitrary energy of the photon (gluon) 
is as yet lacking. In the present paper, based on the technique of Ref. ||, we develop a 
quantum theory of the LPM effect for the whole range of photon (gluon) energies. As 
usual, we treat the medium as a system of uncorrelated static scatterers (atoms). 

In order to set the background we start with the LPM effect in QED. In vacuum the 
evolution of the wave function without radiative corrections of a relativistic electron with 
longitudinal momentum p z in the variable t — (t + z)/2 at (t — z) =const is described by 
the Hamiltonian H = (p 2 + m 2 e )/2^ e ||, where p is a transverse momentum, m e is the 
electron mass, and /i e = p z . Consequently, the electron wave functions at planes z = Z\ 
and z = z 2 are related by 

V>(P2> z 2) = J dP\K e {p 2 , z 2 \p 1 , Zi)tJ)(pi, zx) , (1) 

where p 1 2 are the transverse coordinates, and K e is the Green's function of the two- 
dimensional Schrodinger equation for a particle with mass \x e , times the phase factor 
exp[—iml(z2 — Zx)/2fi e ]. Eq. ([[]) holds for each helicity state. At high energy spin effects 
in interaction of an electron with an atom vanish, and equation analogous to ([[]) holds for 
propagation of an electron through medium as well. The corresponding propagator reads 

K e (p 2 , z 2 \Px, zx) = J Vpexp^i J dz 

where p = dp/dz, U(p, z) is the potential of the medium. 

The interaction of an electron with the photon field generates the radiative correction, 
5K e , to the propagator (Q). To leading order in a = 1/137, 5K e is generated by sequential 
transitions e — > e'7 — > e. The probability of passage of an electron through the target 
without radiation of the photon can be written in terms of K e and SK e as 

P e = 1 + 2Re J dp 1 dp' 1 dp 2 [(5K e (p 2 , z 2 \p x , zi)K* e (p 2 , z 2 \p[, zx)) - (vac)] , (3) 
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where (...) means averaging over the states of the target, and (vac) denotes the vacuum 
5K e K* term. The subtraction of the vacuum term takes into account the renormalization 
of the electron wave function. Eq. (^|) is written for the target with the density independent 
of p, and the points Z\ and z 2 are assumed to be at large distances before and after the 
target, respectively. The initial electron flux is normalized to unity. Evaluation of P e 
allows to determine the probability of radiation of the photon, P 7 , through the unitarity 
relation P e + P 1 — 1. 

The contribution of transitions e — > e'j —>■ e to 5K e can be written as 

1 2 2 2 2 

5K e (p 2 , z 2 \pn zi) = - J dx J d£i J d£ 2 J dT 1 dT 2 g(£i,£,2,x) 

o *i a 

xK e (p 2 , z 2 \t 2 , &)K el (r 2 , 6|ti, £l)K 7 (t 2 , z 2 \t u £i)iT e (Ti, £i|p 1; z x ) . (4) 

Here the indices e' and 7 label the electron and photon propagators for the intermediate 
e'7 state. The transverse masses which enter K e > and K 1 are given by /v = (1 — x)/i e and 
fi-y = xfi e , where x is the light-cone Sudakov variable of the photon. The vertex function 
g(£.i,&,x) equals 

<7(&,£ 2 ,z) = A n/ (x)[v 7 (£ 2 ) - v e ,(£ 2 )] ■ [v 7 (d) - v e ,(ei)] + A s/ (x) , (5) 

where A n f(x) = a [4 — 4x + 2x 2 ]/4x, A s /(x) = am 2 e x/2ii 2 el) v 7 and v e / are the transverse 
velocity operators, which act on the corresponding propagators in Eq. @. The two terms 
in (H) correspond to the e — > e'7 transitions conserving (nf) and changing (sf) the electron 
helicity. 

Making use of Eqs. (0)-(fD we can write the rate of the bremsstrahlung in the following 
differential form (we suppress the vacuum term, which will be recovered in the final 
formula (|15|)) 

dP 22 22 

— - = 2Re J d^i J d^ 2 J dp^p^dridr'^r^T^dp^^i, £ 2 , x) 
^S(p 2 ,p 2} z 2 \t 2 ,t' 2 ^ 2 )M{t 2 ,t' 2 ^ 2 \t 1} t\^ 1 )S(t u t' 1 ^ 1 \p 1 ,p' 1 ,z 1 ) , (6) 
where S and M are defined as 



S(p 2 ,pU2|Pi,PUi) = / T>p e Vp' e exp J d£(p; 



J P'e 2 ) 



®({Pe},{p' e }) 



M{p 2 , p' 2 , £ 2 |pl, P'l, £l) = J VpeVp^Vp^ exp j~ y cf£(/VPe + /l 7 p* - Hepl) 



(7) 



$({p e '},{p e }), (8) 



$({Pi},{P,}) = (expjze / ^[^(0,0-^(0,0]}) 



(9) 



Here // = 2p e (l — x)/m\x is usually called the photon formation length. The boundary 
conditions for trajectories in Eq. (^) are p e (£i, 2 ) = Pi2> Pe(Ci,2) = Pi 2> an d m Eq. (||) 
P e ' 7(^1,2) = Pi 2> P e (£i, 2 ) = Pi 2- Averaging over positions of atoms in Eq. (|9]) yields || 
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$({Pih {P,}) =exp 



din{i)a{\ Pi - Pj \) 
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here n(£) is the target density, and <r(p) is the dipole cross section for scattering of 
the e + e~ pair of transverse size p on a free atom. For the atomic potential 0(r) = 
4n(Za/r) exp(— a/r) (a = IAZ^^tb) cr(p) is given by 

P T S (P" 



a(p) = 87r(ZaaY 



a \a 



For pCa, which will be important in the considered problem, a(p) ~ C(p)p 2 , where 



C(p) = 4vr(Za) 



, /2a\ (l-2 7 ) 

log 7 + 2 



7 = 0.577. 



(12) 



For nuclei of finite radius Ra, Eq. (|1^) holds for p > Ra, and C(p < Ra) = C(Ra) 
Analytical evaluation of the path integral (0) yields || 

Pe \ J 



S(p 2 ,P2,&\Pl,p'l,£,l) 



exp 



\2Ae 



[(Pi - p 2 ) - (pi - p' 2 



den(e)a(|T,(0l) 



/ \2n 
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r s (0 = (Pi-p'i)%# + (p 2 -p 2 )^ T7 ^, A< ±,-t>. 



Introducing in @ the Jacobi variables a = (/vP e ' + P-yP-y)/(pe> + A* 7 ) an d f3 = p e , — p , 
and integrating over the trajectories and p e (£) one can obtain 



M(p 2 ,p2,6|Pi,pi,6 



\ 2 r • 



vPl - P2) ~ (Pi - P2 



2ttA£ 
x J £>/3exp J ijd£ 



zA£' 



f 



/V 7 /3 . ^(Q^(|T m (Q|) 

_ — r t _ 



, (14) 



Tm(U = (Pi - Pi) + (P 2 - P 2 ) Tl— + 



■ (/i £ ' + /i 7 )' 

where /v 7 = /VA^A/V + A*?) = -^e^(l — (-^e is the electron energy). 

Substituting (0), (|14|) into (§), and integrating in (^) over the transverse variables we 
finally obtain (we set —z\ = z 2 = oo, and recover the vacuum term) 



dP oo oo i , > 

-2Re / dfr /rf6exp(-^U(6,6^)^ 7 (0,6|0,6)-^ 7 (0,6|0^ 



-oo ^ 
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V\n=0 
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/v 7 /3 
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n(0a(|/3|x) 
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Thus, we expressed the intensity of the photon radiation through the Green's function of 
the Schrodinger equation with the imaginary potential (|17D . Equation fll5|) is the main 
result of the present paper. 

To proceed with analytical evaluation of the radiation density we take advantage of the 
slow /^-dependence of C(\/3\x) at \f3\x < l/m e which as we will show below are important 
in Eq. (|15|) . Evidently, to logarithmic accuracy we can replace (|17D by the harmonic 
oscillator potential with the frequency 



(1-i) (nC{l3 eff x)x 2 \ 1/2 _ (nC{(3 eff x)x\ l/2 



Here /3 e // is the typical value of \j3\ for trajectories contributing to the radiation density. 
Making use of the oscillator Green's function, after some algebra one can obtain from 
Eq. flT5| ) the intensity of bremsstrahlung per unit length in the infinite medium 



dP 1 
dxdL 



\C{l/m e )_ 



da 

dx 



BH 



Snfiv) + 



da 

dx 



BH 



(19) 



nf \ / sf 

Here rj = lf\£l\, and (da/dx)^ s j- are the Bethe-Heitler cross sections conserving and 
changing the electron helicity. The factors S n f, S s f are given by 



VV2 



dy 



sh y 



exp 



vV2 y 



cos 



+ sin 



(20) 



sin 



(21) 



At small r) S n f ~ 1 — 16r/ 4 /21, and S s f ~ 1 — 31i] A /21. Up to the slowly dependent on 
rj factor C(p e ffx)/C(l/m e ), the suppression of bremsstrahlung at 77 ^> 1 is controlled by 
the asymptotic behavior of the factors (|20D , fl2lf) : S^/ ~ 3/r/\/2, 5*5/ ~ 3tc/2t] 2 . Eqs. (|20|) , 
(PI) allow to estimate /3 e //- The variable of integration in ( P0|) , (21) in terms of A£ in 
Eq. (|15|) equals |A£fl|. Therefore, for typical values of A£ contributing to the integral (|15[) , 
A£ e /j, we have A£ e // ~ min(//, Then, /3 e // can be estimated from the obvious 

Schrodinger diffusion relation /5 e // ~ (2A^ e ////i e / 7 ) 1 / 2 . In the limit of low density, when 
r] — > 0, we get from this relation /3 e // ~ l/m e a;, and the right-hand side of fllTf ) reduces to 
the Bethe-Heitler cross section times the target density. In the soft photon limit (x — > 0) 
at n fixed, Eqs. QTJ)-© yield 



dP 7 
dxdL 



2a 2 Z /2^1og(2a//5 e// x) 



ttE p x 
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where f3, 



eff 



71 



(Za) 2 nE e x 3 log(2/aZ 1 / 3 ) 



-1/4 



. Notice, that (|22| ) differs from the predic- 
tion of Ref. which corresponds to a replacement of j3 e ffX by i?^ in fl2"2"|). 

Let us now briefly consider a gluon bremsstrahlung by a quark interacting with the 
color screened Coulomb centres. Following Refs. |^ we introduce a gluon mass m g = 
1/aQCDi and treat the interaction of a quark with each centre in the Born approximation. 
The derivation of the gluon radiation rate follows closely the QED case. Making use of 



the fact that — T* = Tg (here T q> q are the color generators for quark and antiquark) and of 



the closure over final states of color centres, one can easily show that in QCD one should 
only replace in ( |TT| ) the QED dipole cross section er(|/3|a;) by the QCD three-body cross 
section of scattering of the color singlet qqg system on the color centre, cr 3 0. The latter 
in terms of the color dipole cross section for scattering of a color singlet qq pair, a 2 , is 
given by 

os(l/9|) = I + a 2 (|/3|(l - x))] - la 2 (\(3\x) . (23) 

For the formation length in QCD we have If = 2E q x(l — x)/[m 2 q x 2 + rn? g {l — x)]. The 
oscillator approximation is applicable for evaluation of the radiation density at tjqcd ^ 
1. In this regime (3 e jf <C clqcd, and to logarithmic accuracy the intensity of gluon 
bremsstrahlung can be evaluated making use of the oscillator Green's function with the 
oscillator frequency 

(1-0 f nC 3 ((3 eff ,x) \ 1/2 _ (1-i) ( nC 3 (/3 eff ,x) \ 1/2 
fiQCD " ~K [ to, ) ~^JT\E q x{l-x)) ■ (24) 



Here 



C 3 (f3 eff ,x) = l[C 2 (f3 eff ) +C 2 (/3 e// (l -x))} - \c 2 {(5 eff x) 



and C 2 (p) — (a^nCx/ft) log (2clqcd/p) (Ct is the second order Casimir for the color 
centre) . 

In QCD the limit x — > at the same time corresponds to t]qcd 0. Thus at 
sufficiently small x the Bethe-Heitler regime takes place. On the other hand, at t]qcd 3> 1 
and 1 C 1 the oscillator approximation yields 



dp g _ 9 2 nC T log(2a QCD /Peff ) , , 

dxdL~ as \ 37cE q x3 ' [ } 

with (3 e ff ~ [xa^CTnEq]' 1 / 4 . Eq. (|25| ) disagrees with dP g /dxdL oc x~ 3 predicted in 
Ref. i. 
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